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Abstract. We present a phenomenological study of the phase dynamics of the one-dimensional stochasti¬ 
cally forced Burgers equation, and of the same equation under a Fourier mode reduction on a fractal set. 

We study the connection between coherent structures in real space and the evolution of triads in Fourier 
space. Concerning the one-dimensional case, we find that triad phases show alignments and synchroni¬ 
sations that favour energy fluxes towards small scales -a direct cascade. In addition, strongly dissipative 
real-space structures are associated with entangled correlations amongst the phase precession frequencies 
and the amplitude evolution of Fourier triads. As a result, triad precession frequencies show a non-Gaussian 
distribution with multiple peaks and fat tails, and there is a significant correlation between triad precession 
frequencies and amplitude growth. Links with dynamical systems approach are briefly discussed, such as 
the role of unstable critical points in state space. On the other hand, by reducing the fractal dimension 
D of the underlying Fourier set, we observe: i) a tendency toward a more Gaussian statistics, ii) a loss 
of alignment of triad phases leading to a depletion of the energy flux, and iii) the simultaneous reduction 
of the correlation between the growth of Fourier mode amplitudes and the precession frequencies of triad 
phases. 


1 Introduction 


The physics of extended nonlinear dynamical systems is 
often characterised by fluctuations at different frequencies 
and/or wavelengths. In many cases it is key to disentan¬ 
gle the evolution of amplitudes and phases of such fluc¬ 
tuations [n[2ii3| . In this paper we consider the context of 
nonlinear partial differential equations (PDEs) describing 
the interaction of oscillating Fourier modes (waves) with 
quadratic nonlinearities in general, and the application to 
the one dimensional Burgers equation: 


du 1 du^ d'^u 
dt 2 dx dx'^ 


( 1 ) 


where u(x, t) is the space periodic velocity, is a positive 
parameter (viscosity) and f{x,t) is the external forcing. 
The Burgers equation, see [4] for a review, is a model for 
various nonlinear dissipative systems. It describes a vari¬ 
ety of nonlinear wave phenomena such as acoustic waves 
and plasma physics [S]. When driven by a random forc¬ 
ing it has applications in condensed matter problems like 
interface deposition and growth (see for instance i)- In 
absence of a driving force, the explicit solution provided by 
the Hopf-Cole method wm leads to non-trivial problems 
in the limit of vanishing viscosity, when random initial 


conditions are assumed. This is particularly important for 
the evolution of large-scale structures in the Universe P 
mmmm- Burgers equation was originally conceived 
as a toy model for turbulence and it is frequently used 
as a testing ground for numerical schemes and as a train¬ 
ing ground for developing mathematical tools to study 
Navier-Stokes turbulence and other hydrodynamical or 

Burgers equation represents one of the simplest nonlin¬ 
ear partial differential equations known to display a non¬ 
trivial scaling of the velocity field correlation functions. 
Multiscaling is connected to the tendency to create shocks 
and consequently to increase negative velocity differences 
ArU < 0 and decrease positive ones ArU > 0, where 
ArU = u{x + r,t) — u{x, t). Thus a strongly non-Gaussian 
probability distribution function (PDF) of ArU is observed. 
The presence of strongly localised velocity jumps (shocks) 
in the real space is the fingerprint of Burgers’s dynamics. 
We remark that shocks are the only structures in the flow 
able to dissipate energy in the limit of small viscosity. In 
other words, the energy flux across scales is absorbed only 
by a few strongly localised events in real space. The guid¬ 
ing motivation of this paper is to look for relations between 
the real-space multiscaling properties and the Fourier dy¬ 
namics, a key problem also in Navier-Stokes turbulence 

Ennu. 
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We will work almost exclusively with Fourier mode 
variables Uk{t) rather than the real periodic field u{x,t), 
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where u(x, t) = with x the position in real 

space and k the wavenumbers. Reality of u{x, t) implies 
U-k{t) = ul{t), where * denotes complex conjugation. 
When Galerkin truncations are considered, the range for 
the sum above is reduced: instead of summing over A: S Z 
we sum over a given subset: k G C, where C C Z. Each 
mode is indexed by an integer wavenumber k. The dy¬ 
namical content of mode k is given by its complex valued 
Fourier amplitude, Uk(t). Using this representation, the 
governing PDE can be decomposed into a set of ordinary 
differential equations (ODEs) which describe the individ¬ 
ual evolution of each of the complex amplitudes. To fur¬ 
ther separate the variables of interest an amplitude/phase 
representation will be used: Ufe(t) = where 

o^kit) = \uk{t)\ (amplitude) and ipk{t) = arg[ufc(t)] (phase). 

In Burgers equations the nonlinearity is quadratic so 
interactions appear in triads (groups of 3 modes). The 
key dynamical degrees of freedom consist of the modes’ 
real amplitudes ak{t) along with the triad phases^ i.e. the 
combinations: 

‘A'fe?. fc2 (^) = (^) + '^fc2 (0 - (i) > 

where the wavenumbers ki,k 2 ,k 3 satisfy a ‘closed-triad’ 
condition: ki + k 2 = k^. By simple counting, one obtains 
that many triad phases are not linearly independent; a 
maximal set of linearly independent triad phases is easily 
found once the set of wavenumbers C is known. From here 
on, we will call “state space” the collection of modes’ real 
amplitudes along with a maximal set of linearly indepen¬ 
dent triad phases. 

The motivation for studying the evolution of these key 
dynamical degrees of freedom is that they provide quan¬ 
titative information about the energy exchanges taking 
place in the system. For example, it was shown in |22] for 
the 2D barotropic vorticity equation with periodic bound¬ 
ary conditions that the triad phases j,^{t) not only 
oscillate in time but also precess (or drift). This preces¬ 
sion gives rise to a new frequency, the so-called precession 
frequency, which is approximately equal to the value of 
the zero-mode in time of ^.^(A)- When the precession 
frequency of a given triad is equal to zero or coincides 
with a typical frequency of the system’s oscillation, then 
there is a precession resonance, characterised by strong 
energy transfers. To see how this works, the contribution 
from a given triad to the energy flux towards small scales 
(to be defined below in Section is proportional to the 
time/ensemble average {ak^ak^ciks fc^)) • Two reso¬ 

nant scenarios are possible regarding the phase behaviour, 
i) The triad phase simply oscillates near the value tt/2 
mod 27r. We get a non-zero contribution to the energy 
transfer from this average because ak^Ok^^Ok,^ > 0 by deh- 
nition. This corresponds to precession resonance with zero 
precession frequency: ~ 0. ii) The phase drifts 

with non-zero precession frequency, for example (t) « 

f2t. Thus, the contribution to the flux has the form 

sin(J7t)), which will not vanish if 17 coincides 
with a typical frequency of the time signal of 
This corresponds to precession resonance with non-zero 


precession frequency: ^^) « 17 7 ^ 0. Both resonant 

scenarios can involve multiple triads in synchrony, result¬ 
ing in strong energy transfers across multiple scales of the 
system. This phenomenon leads naturally to the study of 
invariant manifolds within the state space. In |22j . pre¬ 
cession resonances were shown to correspond dynamically 
to trajectories embedded in the unstable manifolds of pe¬ 
riodic orbits. These unstable manifolds have components 
that extend towards large values of amplitudes, ak(t), con¬ 
tributing directly to intense energy-transfer events. 

The aim of this study is twofold. First, we want to 
study the relevance of triad phases and their precession 
frequencies relating to the energy transfers amongst modes, 
focusing on: i) Probability distribution of triad phases, 
joint probability distributions between different triad phases 
and probability distribution of triad phase precession fre¬ 
quencies. ii) Joint probability distributions between pre¬ 
cession frequencies and energy content and also between 
precession frequencies and growth rate of energy trans¬ 
fers. Second, in order to better understand the connections 
between these statistical properties and the presence of 
shock-like structure in the real space, we will study the ef¬ 
fect of the introduction of a quenched disorder given by the 
projection of the dynamics on a fractal Fourier set, where 
preliminary results show a strong depletion of shocks |23j 
with reduced system fractal dimension. We will comple¬ 
ment these studies with dynamical-system interpretations 
in terms of correlations and synchronisation within the 
state space. 


2 Burgers equation in one dimension and less 

Let us start from equation Q written in Fourier space: 

^kiUk-ki - J^k'^Uk + fk: ( 2 ) 

ki 

where energy is dissipated from the viscous term and in¬ 
jected by the external forcing. The core of the dynamics is 
in the quadratic non-linear convolution, a term that glob¬ 
ally conserves energy by redistributing it amongst Fourier 
modes. Let us define the stationary Energy spectrum as: 

{ukUk') = EkS{k + k'), (3) 

where (•) denotes an average over the statistically-steady 
ensemble. The energy flux in the positive direction across 

a wavenumber k is defined as n{k) = —^ 

From ([^ and ([^ it is easy to realise that the contribution 
from the nonlinear interactions to this energy flux is given 
by the “triadic terms”: 

k 00 

n{k)=Y^ ki^{{uk,Uk2U*k,+k2)} ^ ( 4 ) 

ki— 1 k2 — — oo 

where 3 { 0 } denotes the imaginary part of z. It is impor¬ 
tant to point out that the energy exchange produced by 
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Fig. 1. Left: Snapshots of real space velocity profiles as a function of space and fractal Fourier dimension: D = 1 (red), D = 0.99 
(blue), D = 0.97 (grey) and D = 0.95 (black). Right: Mean energy spectra in log-log scale, obtained for the same four fractal 
dimensions D and using the same colour code as in the left panel. The black solid line is the scaling slope k~^ produced by the 
undecimated {D = 1.0) Burgers equation. 


the nonlinear term inside a single triad is conservative. 
Indeed, the sum of the energy variation at three wave- 
numbers with ki k 2 = fca is always zero. This implies 
that the energy flux across a scale k is only given by the 
energy exchanges inside triads which have at least one 
wave number higher than k and one smaller than k. Hence, 
7T(oo) is always equal to zero [24l[^l26l[27] . In particular, 
at those scales where neither the forcing mechanism nor 
the viscous term play any relevant role, the dynamics are 
fully dominated by the non-linear evolution, the so-called 
“inertial-range” of scales. For smaller and smaller viscos¬ 
ity, the solution of Burgers equations develops a velocity 
configuration with sharper and sharper shock-like struc¬ 
tures (see fig. [^, characterised by a power-law spectrum 
E{k) oc k~'^ (see [4]). 


2.1 Phases and fluxes 

Focusing on the non-linear term of the Burgers equation 
(11) and making the following substitution 

Mfc(t) = afc(i)e*'^'“‘'*\ ak{t) = \uk{t)\, (5) 

and recalling the definition of triad phase 

+ 0fc2W -'/'fcaW: (6) 

we get the evolution for the amplitude ak and phase (j)k of 
the surviving modes: 

k - 


= - 0 


0^ki0^k2 


2 CLk 

ki ,/C 2 


Ok^Ok^ cos{ipl^^k2) 5k2+k2,k- ( 8 ) 


Note that the factor 9k in the latter equations is the 
projector in front of the nonlinear term in eq. 0 , while 


9ki and Ok^ represent the fact that also the other two 
modes satisfying the triadic condition are decimated. 

It is evident from eqs. (j^, (j^ and (j^ that ^kl k^^ 
along with the Uk, satisfy a closed system of equations 
which determines the nonlinear dynamics. In particular, 
the individual phases (pk can be obtained by quadratures 
from eq. (|^. 

Applying substitutions ([^~® to the velocity field in 
the energy flux equation, eq. (BVwe get: 

k 00 

n{k)^Y^ (ak.ak^ak.+k^siniipll+^'l^)'^ 

ki—1 k2 — — oo 

X ki0kiOk2^ki-\-k2- (^) 

As expected, the energy flux depends on the phases via 
ocily. The crucial observation is that the sign of the 
contributions to this energy flux depends only on the val¬ 
ues of the triad phases, because all other factors appearing 
in ([^ are positive, including the real amplitudes ak- Since 
the sign of the energy flux determines the cascade direc¬ 
tion, it is very important to analyse these contributions in 
detail. 

Define fca by the triad relation fci + ^2 = ^3 ■ Following 
an analogous method to the detailed energy balance of tri¬ 
ads we deduce that non-zero contributions 

to the energy flux ([^ come from triads that have at least 
1 mode within the interval [—fc, fc] and at least 1 mode 
outside it. Combining these contributions we obtain 

k 00 

n{k)=Y 2/ci (|afciafc,afc3sin(v?*3 , 

fci—1fca—fc+l 

where we have used the reality of the original field, which 
implies 4>-kj = ~4’kj and a-kj = o-kj- We remark the 
wavenumber ordering: 0 < fci < fca, 0 < k 2 < k^ and 
^3 = ^1 + k2- 
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In summary, the sign of each individual triad’s con¬ 
tribution depends solely on the sign of sin(i^^|J where 
ki,k 2 , fca is the “ordered” version of that triad, with = 
fci -I- k 2 and kj > 0. As a result, when is close to 

I -I- 2mr, n G Z, the direct-cascade e nerg y flux is max¬ 
imised (at fixed amplitudes). In Section [3.1| we present sta¬ 
tistical analyses that show that triad phases tend to clus¬ 
ter near these values, involving in addition a synchronisa¬ 
tion amongst different triads, so that several terms con¬ 
tribute with a positive sign, leading to direct-cascade flux 
whose intensity depends on the dimension of the Fourier 
fractal set where modes live. 


Table 1. D: dimension of the Fourier Fractal set where modes 
live. N: number of collocation points. %{D): percentage of dec¬ 
imated wave numbers, where the first value is related to the 
lower resolution used while the second value is related to the 
higher resolution, v: value of the kinematic viscosity, kf. forced 
wavenumbers. St: Numerical Time step used in the temporal 
evolution. 


D 

N 

%(D) 

V 

k 

/ 

St 

1 

216 _ 219 

0 

m 

1 

O 

i-H 

00 

[1 

8 ] 

5.5- 10“® 

0.99 

216 _ 219 

8-10 

2.5 • 10“® 

[1 

8 ] 

2.3- 10“® 

0.97 

216 _ 219 

23-27 

9-10“® 

[1 

8 ] 

2 .0- 10“® 

0.95 

216 _ 219 

36-40 

5-10“® 

[1 

8 ] 

1.7- 10“® 


2.2 Burgers equation on fractal Fourier sets 


In order to study the effects of the quenched noise on the 
triadic interactions in Fourier space we use a projection 
method to fractal sets. This projection method was orig¬ 
inally introduced in [55] to study the inverse energy cas¬ 
cade in 2D Navier-Stokes equations (see |53| for a review). 
Later, it was exploited to study the statistical properties of 
small-scale fluctuations in the 3D Navier-Stokes equations 
[30] and ID Burgers equation [23] ■ In these two papers, a 
tendency towards more regular Gaussian statistics for the 
small-scale velocity field was observed at lower values of 
fractal dimension D, indicating that the disorder leads 
to an energy transfer regime characterised by self-similar 
fluctuations. 

The projection on a fractal Fourier set is done using a 
decimation operator Pd acting on the field u{x, t) as 


v{x,t) = PDu{x,t) = (10) 

k 


Here 0/j, fc S N are independently chosen random numbers 
such that = 1, with probability and Ok = 0, with 
probability l — hk- We also impose the condition that 9k = 
6_k to ensure the reality of the field at all times. Choosing 
hk = k^~^, with 0 < D < 1, we introduce a quenched 
disorder which randomly suppresses modes on the Fourier 
lattice and ensures that on average we have N{k) oc k^ 
surviving modes inside an interval of length k around the 
origin. Applying this projector to the Burgers equation, 
we can then write the decimated PDF as 


dv 1 dv^ d'^v 

m ^ P°l)x ^ 


+ PDf- 


( 11 ) 


We performed a set of numerical simulations of eq. (Ill 
by changing the dimension between D — 1 and D = 0.95. 
We chose the forcing to be Gaussian and white-in-time 


(/>i,D)/(fc2,t2)) =2 /o|fc/|-i<5(ti-t2)-5(fci + fc 2 ), (12) 


acting only at large scales, the forced wave-numbers kf are 
all the non-decimated modes in the range [1:8]. For the 
numerical simulation we used Adams-Bashforth schemes 
of fourth order, with a number of collocation points ranged 
between N = 2^® and N = 2^®, and a time step 6t ~ 10“®. 


See Table[2for more details about the numerical data. Fig¬ 
ure]^ shows a snapshot of the velocity field for four differ¬ 
ent values of fractal dimension D, together with the mean 
energy spectra. It is important to mention that in fig. [^ 
the spectra have no gaps because we have further per¬ 
formed an average over different quenched fractal masks 
[23] . As one can see, shocks are present for all fractal di¬ 
mensions, while the smooth (for D = 1) ramps connect¬ 
ing two shocks develop small-scale fluctuations that be¬ 
come more and more pronounced as the fractal dimension 
D decreases. It is important to remark that even if the 
disorder produces a non-trivial change in the spectrum, a 
power-law behaviour is always present |23| . Similar results 
have also been observed in the three dimensional Navier- 
Stokes equations [50] ■ A scaling of k~'^ is recovered in the 
undecimated D = 1 case. 


3 Numerical results 

3.1 Phase analysis 

We present a phenomenological study of Burgers triad 
phase dynamics. Let us first consider the case when we 
define all phases in the same periodic range, (t) G 

[—7r,7r], i.e., we omit the phase precessions for the mo¬ 
ment. In hg. 0 (left panel) we show the PDF of 
for three difteent triads composed of wave numbers be¬ 
longing to the inertial range and ordered in the usual way 
0 < fci < /c2 < ^3 with k^ = ki k 2 '. 

Ti^ = [100;150;250], 

T2^ [/ci;A:2;fc3] = [200;250;450], 

Tls ^ [/ci; fe; ks] = [300; 350; 650], 

at D = 1 and D = 0.95. The undecimated case D = 1 
shows a more marked preference for triad phase align¬ 
ments about 7r/2, which maximises the energy flux be¬ 
cause flux contributions are proportional to sin((/?^’[ f,^) 
(see eq. (©)■ The PDFs are peaked at 7r/2 and have a 
minimum at 37r/2, which suggests that the energy flux is, 
on average, directed from large to small scales. This would 
be true provided one assumes that real amplitudes Ofc take 
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Fig. 2. Left: Histograms of triad phase (mod 2?! with offset —tt) computed during the temporal evolution of different triads 
in the inertial range. Ti : [k\-,k 2 ',ki\ = [100; 150; 250] (black solid line), T 2 : [ki-,k 2 ',k^\ = [200; 250; 450] (red dashed line), and 
T 3 : ]fei;fc 2 ;fc 3 ] = [300; 350; 650] (blue dotted line). To aid comparison, the undecimated case {D = 1) and decimated case 
[D = 0.95) are plotted together. Right: Same type of histograms but over all 160801 triads composed by wave numbers inside 
the range 100 < k < 1000, and for the same set of values of fractal Fourier dimension D and colour code as in fig. 


higher or similar values when phases are near 7r/2 than 
when phases are near 37 r/ 2 . This assumption is confirmed 
by simulations, so we have a direct cascade of energy as in 
3D Navier-Stokes turbulence. On the other hand, for the 
case D = 0.95 we no longer have a marked peak at 7 r/ 2 , 
especially when we consider the triad T 3 which is the one 
containing higher wavenumbers. This is the first evidence 
that the introduction of the quenched disorder strongly 
perturbs the coherence of the energy transfer mechanism, 
leading to a much slower direct cascade. In hg. (right 
panel) we show a similar study but this time involving 
over 160 thousand triads with wavenumbers in the range 
[100,1000]. The features just discussed for the selected 
triads are again there, which shows that most triads are 
behaving in basically the same way. 

A qualitative view of the correlation amongst different 
triad phases throughout the system’s evolution is pro¬ 
posed in fig. Here a scatter plot between the phase 
of triads Ti and T 2 is shown, where the phase angle is 
not constrained to be on a periodic domain. From the left 
panel (undecimated Burgers equations D = 1) it is clear 
that a correlation is present between the two phases. A 
lattice-like distribution with the presence of many clus¬ 
ters indicates that for long periods the two phases are 
oscillating around a fixed value before jumping simulta¬ 
neously toward other values. As expected, these clusters 
exist around values of 7 r/ 2 -|- 2 n 7 r, n G Z. This synchroniza¬ 
tion is directly linked to the results shown in Hg. where 
we saw a peak in the PDF for phase values at 7 r/ 2 . Since 
most triad phases spend most of their time near 7r/2, it 
is natural that several triads will be simultaneously near 
those values for some of the time. In contrast, results for 
D — 0.99 and 0.95 in fig.[^ centre and right panels, show 
that the synchronization amongst different triad phases is 
washed out when the quenched disorder is introduced. 

Remarkably, synchronisation of triad phases has a di¬ 
rect effect on the energy flux. The more triads are syn¬ 


chronised near 7r/2 -|- 2n7r, the more positive terms add 
together to contribute to the flux in eq. (1^, leading to 
a large flux. Thus, loss of synchronisation diminishes the 
flux via cancellations of positive and negative terms. 

3.2 Precession analysis 

The “drift” or “precession” of triad phases made evident 
from fig. 1 ^ suggests that we study the distribution of pre¬ 
cession frequencies of the various triad phases in the sys¬ 
tem, as this distribution can have an impact on the energy 
fluxes via precession resonance [22] . In order to analyse the 
phase-precession statistical and dynamical properties it is 
important to Alter out the high-frequency fluctuations in¬ 
troduced by the evolution of each of the phase variables 
contributing to the triad phase. Therefore, let us define 
the triad phase precession averaged on a time window At 
as 

.. t+At 

H ‘^klk2 (^0 > (13) 

ti=t 

where ti denote data points from the numerical integra¬ 
tion and the time window. At = N5t^ where 5t is the 
numerical time step and is a large enough integer to 
remove the noise at the dissipative scales. But At should 
be chosen small enough not to kill all temporal correla¬ 
tions in the system. In practice we choose At of the order 
of the Kolmogorov time microscale trj = « 10“^. 

The instantaneous phase precession, is defined 

in terms of the finite-difference increments of individual 
phase derivatives ipk(t): 

Ukjtj + St) 

UkiU) 

(14) 

The probability density functions (PDFs) of {t, At) 
are produced by considering all 160801 triads constructed 


^k(ii) — 


<Pk[ti + St) — (j)k{ti) 
St 
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Fig. 3. Parametric plots of the phase evolution for two different triads in the inertial range: Ti : [fei; ^ 2 ; k^] = [100; 150; 250] and 
T 2 : [fci; k 2 ’, ks] = [200; 250; 450], showing synchronisation events near ^ + 2mT , n € Z. Left: Fractal Fourier dimension D = 1. 
Grid lines denote triad phase values at ^ + 2nn. Centre: Fractal Fourier dimension D = 0.99. Right: Fractal Fourier dimension 
D = 0.95. 


from modes inside an interval in the inertial range: k € 
Ik = [100,1000]. In fig.|^we show the corresponding PDFs 
for three different time windows: At = 1000 = 0.51,, 

(top row); At = 2000= 1.01,, (middle row); At = 
6000= Z.Qtri (bottom row). Furthermore, we also com¬ 
pare the results obtained from the original undecimated 
{D = 1) Burgers equation (left), with two different deci¬ 
mated systems: D = 0.99 (centre) and D — 0.95 (right). 
We first note the appearance of secondary peaks in the 
PDFs, corresponding to non-zero precession frequency val¬ 
ues. The number of peaks increases as the time window 
At increases because the triad phase has had enough time 
to make several jumps between adjacent trapping regions, 
namely from 7r/2 -|- 2n7r to 7r/2 -|- 2n'7r, with jn — n'\ = 1, 
in apparently random sequences (figure not shown, but 
the jumpy behaviour is evident from the scatter plots of 
fig. |3]). Similarly, the peaks become more pronounced as 
At increases because the longer averaging window leads 
to a suppression of oscillations in around the trap¬ 

ping regions 7r/2 -|- 2n7r. This explanation also applies to 
the secondary peaks where the time window At includes 
both the oscillations of referred to above but also a 

number of jumps in towards other trapping regions. 

As At is increased the contributions of the oscillations are 
suppressed but the contributions from the jumps will accu¬ 
mulate. The second interesting feature is the sudden dis¬ 
appearance of the secondary precession-frequency peaks 
as we decimate the system. At D = 0.99 they are barely 
detectable and at D = 0.95 they have been wiped out. 
This is in line with the results presented in figs. and 
showing a transition from a structured phase evolution for 
D = 1 to a disordered distribution when decimation is ap¬ 
plied. 

In order to explore the energy flux it is also important to 
examine the correlation between evolution of triad preces¬ 
sion and triad amplitude. For the same subset of triads, 
ki + k 2 = fca in the range Ik, we computed two different 


joint PDFs, corresponding to the correlation amongst the 
triad precession I 2 klk 2 temporal average, 

over At, of: i) the modulus of the product of the mode 
amplitudes, 

, (15) 

and ii) the modulus of the product of the logarithmic time 
derivatives of the mode amplitudes, 

^k2 ^^3 

(kki Ojk2 ^^3 

The motivation for choosing these quantities is that they 
provide a fair comparison weight, amongst different triads, 
with respect to energy content i) and flux of energy ii). 

In fig. [^we show the resulting joint PDFs for a time 
window At = 3.0together with its dependency on the 
fractal dimension D. From these PDFs one can again de¬ 
tect the presence of secondary correlation peaks in the 
undecimated case only {D = 1). The introduction of the 
quenched disorder destroys the cross-correlations amongst 
amplitudes and phase precessions. 

A closer inspection of fig. reveals a connection with 
a dynamical systems approach, which deserves some men¬ 
tion here. By definition, a critical point of system ([^ is 
a point in the state space where all time derivatives dk 
and <fklk 2 equal to zero simultaneously, for all modes 
and triads. The most common type of critical point is 
the so-called “unstable”: the state of the system may ap¬ 
proach such a point along some directions, getting close 
to it but eventually separating from it along other direc¬ 
tions. Thus, active energy exchanges are expected if the 
state of the system gets close to an unstable critical point. 
The evidence from fig. is that the state of the Burgers 
system visits unstable critical points regularly during the 
evolution: notice the tail near dk = 0, <fklk 2 =0 3-1 1^6 
lower part of the bottom left panel. Moreover, this tail is 
correlated, by definition, with the tail at the upper part 
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Fig. 4. Probability density functions (PDFs) of the time-averaged triad precessions At) (over time window At, see 

eq. (13l). The PDFs are computed over all 160,801 triads composed by wave numbers inside the range 100 < k < 1000 and 
for the temporal evolution 0 < t < 10® t,,. Top row: At = 100051 = 0.51,,. Middle row: Zil = 200051 = 1.01,,. Bottom row: 
zil = 600051 = 3.01,,. The three panels in each row represent different fractal dimensions. Left: D = 1. Centre: D = 0.99. Right: 
D = 0.95. 


of the top left panel, corresponding to events of maxi¬ 
mum amplitude, i.e., high energy transfers. Finally, fig. 
suggests that these critical points are distributed prefer¬ 
entially near values of the triad phases given by = 

7r/2-|-2n7r. In contrast, in the fractal decimated case {D = 
0.95; see fig.[^ right panels) the above-mentioned tails lose 
coherence, which could mean that the state visits critical 
points less often or that less critical points exist. These 


and other implications will be studied in detail in subse¬ 
quent work. 


4 Conclusions 

In this work we have presented numerical results about 
the evolution of the phases of the triadic structures re¬ 
sponsible for the energy transport in the Fourier repre¬ 
sentation of ID Burgers equation. Figurej^shows that the 






















































































































































Michele Buzzicotti et al.: Phase and precession evolution in the Burgers equation 



Fig. 5. Joint PDFs for time window At = 3.01,,, for the same set of triads as in fig. [d] Top row: PDFs as functions of triad 
precession frequency and product of triad mode amplitudes. Bottom row: PDFs as functions of triad precession frequency and 
product of logarithmic time derivative of triad mode amplitudes. Left: D = 1. Right: D — 0.95. 


PDFs of these phases are characterised by marked peaks 
at 7r/2+2n7r, consistent with the presence of a well-defined 
energy flux directed towards the small scales. In addition, 
strong correlations between the phase evolution of differ¬ 
ent triads are presented in fig. These correlations hint 
to the entangled interactions and synchronisation devel¬ 
oped by the Burgers nonlinearity amongst all the Fourier 
modes in order to provide the well known energy cascade. 

The PDF of precession frequencies, fig. shows a 
structure of secondary peaks with significant probabil¬ 
ity, circumscribed by a non-Gaussian (fat-tailed) envelope. 
This non-trivial feature motivates the search for higher- 
order precession resonances, responsible for perturbations 
in the energy flux across scales. These resonances occur 
when the precession frequency ^ triad phase co¬ 

incides with a typical nonlinear frequency of the triad’s 
mode amplitudes |22j . Such a search requires the intro¬ 
duction of a tuning parameter. Work in this direction will 
be reported elsewhere. 


Another important aspect of this paper consists of 
the study of the effects produced by the introduction of 
fractal-Fourier decimation for the ID Burgers equation. 
We found clear evidence that decimation drives the sys¬ 
tem to a loss of temporal and spatial correlations among 
the surviving triad amplitudes and phases, as figs. [^I^and 
[^show. In real space, as the dimension D is lowered, the 
velocity held tends to develop huctuations in the ramps 
between shocks. These increasing huctuations could be 
linked to the loss of correlation/synchronisation amongst 
the triad amplitudes and phases. In particular, the signif¬ 
icant correlation between triad precession frequencies and 
amplitude growth in the undecimated case (D = 1) is lost 
as the dimension D is lowered (see hg. [^. 

Looking forward, the dynamical systems point of view 
is worthy of further detailed investigation. We need to un¬ 
derstand how phases organise and synchronise during the 
dynamical evolution, what is the distribution of critical 
points in the state space, and how these features change 
when we reduce the fractal dimension D. In particular, the 
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alternation of trapping events and jumping events experi¬ 
enced by the triad phases in fig. suggests the possibility 
to develop a model for phase evolution based on a biased 
random-walk with waiting times |31j , which in turn could 
explain the fat tails seen in fig. Finally, we stress that 
very similar studies can be carried on for the much more 
complicated case of turbulence in the three dimensional 
Navier-Stokes equations. Work in this direction will be re¬ 
ported elsewhere. 
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